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We calculate the heat transfer between electrons to acoustic and optical phonons in monolayer 
and bilayer graphene (MLG and BLG) within the quasiequilibrium approximation. For acoustic 
phonons, we show how the temperature-power laws of the electron-phonon heat current for BLG 
differ from those previously derived for MLG and note that the high-temperature (neutral-regime) 
power laws for MLG and BLG are also different, with a weaker dependence on the electronic tem- 
perature in the latter. In the general case we evaluate the heat current numerically. We suggest 
that a measurement of the heat current could be used for an experimental determination of the 
electron-acoustic phonon coupling constants, which are not accurately known. However, in a typical 
experiment heat dissipation by electrons at very low temperatures is dominated by diffusion, and 
we estimate the crossover temperature at which acoustic-phonon coupling takes over in a sample 
with Joule heating. At even higher temperatures optical phonons begin to dominate. We study 
some examples of potentially relevant types of optical modes, including in particular the intrinsic 
in-plane modes, and additionally the remote surface phonons of a possible dielectric substrate. 

PACS numbers: 73.22.Pr,73.63.-b,72.20.Ht 



I. INTRODUCTION 

In the presence of Joule heating, the static temper- 
ature of the conduction electrons of a metallic system 
is set by the heat balance between the heating and en- 
ergy relaxation. The latter is caused either by electron 
diffusion away from the heated region or by the energy 
transfer to the lattice via the electron-phonon coupling. 
If the lattice can transfer its energy effectively enough 
to the underlying substrate, a hot-electron situation may 
be reached, where the electronic temperature T e is con- 
siderably higher than that of the (acoustic) phonons, 
T oc P2lMore generally, such a situation with well defined 
separate electron and phonon temperatures is known as 
quasiequilibrium. In analyzing experiments or devices 
that involve heating (or cooling's) effects of this kind, it 
is essential to be able to model the heat currents between 
the various subsystems!^ This is particularly important 
in two-dimensional graphene, which has a small thermal 
volume and is thus easily overheated. Of central interest 
here is the heat current between the electrons and the 
phonons, on which we concentrate in this paper. 

In metals at sufficiently low temperatures only acoustic 
phonons are relevant, and of those the longitudinal (LA) 
modes have the strongest coupling constants. At very low 
temperature the power transferred between the electrons 
and LA phonons (assuming quasiequilibrium) is typically 
of the form 



(1) 



where Vd is the (/-dimensional volume and £ is a coupling 
constant. This form is well known in the case of simple 
three-dimensional (3-D) metals, where the exponent is 
6 = 5 (Rets. |2"I5|) . In lower-dimensional and disordered 
systems different exponents have been f ound . In disor- 
dered thin films 5 = 6 has been observedp^ and in thin 
doubly clamped metallic beam, where the electrons re- 



main 3-D but vibrations are one-dimensional, (5 = 3 and 
Vd is the length of the beamP 

The electron cooling power and temperature relaxation 
times were recently studied theoretically also for mono- 
layer graphene (MLG) pHUD Graphene is a single sheet of 
graphite, i.e. & two-dimensional (2-D) honeycomb lat- 
tice of carbon, which has (semi)metallic properties! 11 * 12 ! 
In its monolayer form the valence and conduction bands 
touch at the two K points of the Brillouin zone (or Dirac 
points), around which the electron dispersion relation is 
conical. In bilayer graphene (BLG) the dispersion is ap- 
proximately parabolic!^ The most important difference 
of graphene to other metals or semiconductors is the "chi- 
ral" character of the charge carriers. However, close to 
the charge neutrality point the properties of graphene 
differ from other metals already by the special forms of 
the dispersion relations. In Ref. [5] it was found that at 
low temperature the power transferred to in-plane LA 
phonons in MLG is of the usual form ([!]) with S = 4 [see 
Eq. (15) below]. In Ref. [9j on the other hand, it was 



shown that in the neutral (or high-temperature) regime 
of MLG the power has a form asymmetric in T e and T ac , 



V d g(n,T e )(T e - T ac ), 



(2) 



where Vd = A, the area of the sample, g(fi, T e ) is a func- 
tion which we specify in Eq. (17) below, and \i the chem- 



ical potential measured from the Dirac point. 

In this paper we revisit the problem of electron-phonon 
power transfer in MLG and consider also the case of BLG 
in the parabolic-band approximation. We show that the 
conditions for the validity of the "low-temperature" and 
"high-temperature" results for MLG mentioned above 
are roughly T e , T ac < Tbg.mlg and T e , T ac > T bg ,mlg, 
respectively, where Tbg.mlg = 2(c/«)|ju|/fes is the 
Bloch-Griineisen temperature of MLG. Here »«lx 10 6 
m/s is the Fermi speed in MLG and c ps 2 x 10 4 m/s 
<C v is the speed of sound. These two limits have 
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been termed the Bloch-Griineisen (BG) and equipar- 
tition (EP) regimes, respectively.^ In these limits we 
recover the previously derived results [Eqs. p"5| ) and 
(17)], but in addition to studying the limits analytically, 



we evaluate the power numerically in the crossover re- 
gion. We then repe at the same analysis for BLG, where 
7bg,blg = 2 (c/«)V7iHAb, with -f X the interlayer cou- 
pling amplitude. Also for BLG, in the BG regime the 
usual power-law form ([!]) for acoustic phonons is found, 
with <5 = 4, while the result in the opposite regime has 
the form |2]) [Eqs. ( fl9| ) and (21)]. In fact the exponent 
8 = 4 is independent of the chirality of the carriers, since 
it appears that in non-disordered systems at low temper- 
atures quite generally 8 = d + 2, where d is the smaller of 
the dimensions of the electron and the phonon systems. 
However, the coupling constants E for MLG and BLG de- 
pend rather differently on doping, with E oc \fi\ in MLG, 
and E oc l/-J\jT\ in BLG. These differences in E can 
be understood based on the electron-phonon relaxation 
rates and electronic specific heats. In the EP regime 
we find that the temperature dependence of g((i, T e ) is 
weaker in BLG {g ~ T e ) than in MLG (g ~ T e 4 ). 

At very low temperatures the most efficient energy re- 
laxation mechanism for electrons is diffusion. In order to 
measure Q e ~ac f° r determining the coupling constant E, 
the temperature should be higher than a certain value 



which we estimate in Sec. VII for a typical experimen- 
tal situation. On the other hand, at even higher tem- 
peratures the optical phonon power becomes dominant. 
To estimate also the temperature for this crossover, we 
consider some simplified models for optical phonons that 
may be relevant. The most obvious ones are the intrin- 
sic in-plane longitudinal (LO) or transverse (TO) modes. 
For graphene on a dielectric substrate, the surface opti- 
cal phonons of the dielectric must also be considered! 15 * 16 ! 
We find that the latter can begin to dominate the energy 
relaxation already at much lower temperatures than the 
intrinsic phonons. 

In suspended graphene, in addition to the LA phonons, 
also out-of-plane acoustic (flexural) modes should be 
considered P21 However, these are disregarded here as 
their coupling to electrons is only of second order in the 
displacements and since they should not be important 
for graphene on a substrate. On the other hand, we also 
do not discuss explicitly the coupling of the graphene 
phonons to those of the substrate, or the heat trans- 
fer (Kapitza resistance) between them, 18 which, although 
an important part of the heat-balance problem, depends 
on the details of the interface and is difficult to model 
microscopically. This coupling may to some extent af- 
fect the acoustic phonon dispersions, and the relevant 
LA phonons should possibly be understood as collective 
surface acoustic modes of the coupled graphene-substrate 
system. Finally, we do not take into account effects of 
impurities, or discuss explicitly electron-electron interac- 
tions, which are assumed to be strong enough to keep the 
electrons in quasiequilibrium. 

This paper is organized as follows. In Sec. In] we sketch 



our Boltzmann-Golden Rule approach for calculating the 
heat current and briefly discuss the quasiequilibrium ap- 
proximation. Next we give a brief reminder on the low- 
energy electronic structure of graphene in Sec. |III| and 
then continue to consider the coupling of the electrons 
to LA phonons via a deformation-potential approach in 
Secs. |IV| and|V] In Sec. |VI| we describe our models for the 
optical phonons and their coupling to electrons, and the 
calculation of the corresponding heat currents. Finally, 
we discuss the experimental consequences of our results 
and the crossover temperatures in Sec. |VII| 



II. ELECTRON-PHONON HEAT TRANSFER 
FROM THE BOLTZMANN THEORY AND 
FERMI GOLDEN RULE 

Below we derive a general expression for the electron- 
phonon heat current by employing a Boltzmann collision 
integral in the Fermi Golden Rule approximation. In 
what follows, we disregard the spin and valley indices, 
which only appear as an additional degeneracy factor 
g e = 4 in the results. 

The position-independent Boltzmann equation de- 
scribing the occupation probability f% of the electron 
excitation with momentum hk in band a is 



dtfk — Se-ph(fk)- 
Here the collision integral is given by 



(3) 



P/3 (4) 



where the Golden-Rule scattering rates are 



Wka^p/3 
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[(nl + l)8 k , p+q 8(ef-^ q ) 



(5) 



In these ri^ is the distribution function of phonons with 
momentum hq and band index 7, and are the elec- 
tron and phonon excitation energies, and we defined 
= £fc — £p- We assume the electron-phonon coupling 
constants to satisfy the symmetries W^t'Z = w pk'q = 
w kp -q- The power with which the phonons cool the 
electrons (i.e. the electron-phonon heat current) is de- 
fined by 

Q = -9t^2^fk = -^2 e k S e- P h(fk)- (6) 



ka. 



ka- 



It is convenient to divide the power into two terms ac- 
cording to whether they describe induced (stimulated) or 
spontaneous processes.^ Thus we find 



— Q'md ~T" Q ; 



spont j 



(7) 



where 

Qind 



; EEE 

q7 p/3 fca 



spont 



a0 



(8a) 



q7 p/3 fca 



X 71 



■e(^)Sk,p+ q S(e^-uj^). 



(8b) 

Here and in the following we assume that the system is 
in quasiequilibrium where electrons are described by the 
electron temperature T e and chemical potential and 
phonons by the temperature T p h- Thus ffj! = f{e%) and 
nl = n ph {tol) where f(E) = {exp[(E - ^)/k B T e ] + l}" 1 
is the Fermi function and Tii(w) = [exp(w / kTi) — 
the Bose function at temperature Tj. To arrive at the 
form m we have applied /(e£)[l - /(eg)] = - [/(eg) - 



t fcp/ 1 



which is valid in equilibrium. 



The only difference between Qspont and Qi n d is there- 
fore the sign and the temperature that enters the Bose 
function. Writing Q ind = Qmd(faT e ,T p h), we have 



Q spont (Mj -^e) — Qind (M; T e : T e ) • 



(9) 



Due to this symmetry, it is only necessary to calculate 
Qind explicitly. It is easy to see that the net power trans- 
fer is zero if T e =T p y l . 

Besides the heat current Q, another useful quantity is 
the linear-response thermal conductance, defined as 



d 



dAT 



Q[H,T- 



AT 



AT 
~2~ 



> (io) 

AT=0 

T=(T e + 



where Q = Q(fi, T e , T ac ) and AT = T e - T ac , 
T ac )/2. This is more convenient for characterizing the 
strength of the electron-phonon coupling, because it only 
depends on a single temperature. 

A more physical reason for concentrating on the linear- 
response regime comes from the assumption of quasiequi- 
librium. The latter is a commonly used approximation 
that avoids the need to include the electron-electron col- 
lision integral and an equally detailed treatment of the 
phonon kineticsP Electrons remain in internal equilib- 
rium if the electron-electron interactions are sufficiently 
strong.^ (For an example of a case where this is not 
valid, see Ref. f20l) The same applies to phonons if their 
relaxation is strong enough, which is typically the case for 
systems on a substrate. However, although these condi- 
tions might not always be satisfied far from equilibrium, 
we do not expect them to be a serious concern in the 
linear response regime. 



III. ELECTRONIC STRUCTURE OF 
MONOLAYER AND BILAYER GRAPHENE 

Here we give a very brief reminder on the low- 
energy band structure of graphene. Figure [l] shows 




FIG. 1: (Color online) Geometrical structure of bilayer 
graphene with AB, or Bernal, stacking and the Slonczewski- 
Weiss-McClure band parameters 70,1, 3, 4. The Al, Bl, A2, 
and B2 atoms of the two sublattices and layers are indicated. 
The nearest-neighbor distance is a ~ 0.14 nm, and the inter- 
layer distance is d « 0.37 nm. 



the geometrical structure of BLG. For a tight-binding 
description we emp loy the Slonczewski- Weiss- McClure 
parameters! 12 l 21 ti s5T j-^ w hat follows we only use the pa- 
rameters 70 « 3 eV, 71 « 0.4 eV and set 73 = 74 = 
0, which is a reasonable approximation at not-too-low 
energies.^ MLG is obtained by setting also 71 = and 
concentrating on one of the layers. The unit cell then 
has two atoms and the Hamiltonian for each of the K 
point of the reciprocal space is a 2 x 2 matrix, whose 
eigenvalues give the conical "Dirac" spectrum. In the 
case of BLG the unit cell has four atoms, and thus the 
low-energy Hamiltonian is a 4 x 4 matrix. However, us- 
ing perturbation theory to leading order in I/71, it may 
be replaced with an effective 2x2 Hamilt onian in the 
basis of the "uncoupled" atoms A2 and Bl! 12 ' 13 ' 24 ' The 
resulting dispersion relation is parabolic. 

Using this parabolic-band approximation for BLG, the 
dispersion relations of the valence (v) and conductance 
(c) bands and the corresponding eigenstates may be writ- 
ten collectively as 



9a(&«*)77l' 



V2 



1, =bs a e 



■ n id 



(11) 



where n = 1 for MLG, and n — 2 for BLG, s a — ±1 for 
a = c,v, 4>f. — arctan(fc J/ / k x ) , k = \k\, and Hv = Sajo/2, 
with v rs 10 6 m/s. The upper and lower signs correspond 
to the two different K points (valleys), relative to which 
the wave vectors k are counted. The zero of energy is set 
at the charge neutrality point. 

These low-energy results are sufficient for the dis- 
cussion of interactions of the electrons with acous- 
tic phonons. Atomic-scale scatterers such as short- 
wavelength optical phonons can also couple the two K 
points J^^SI font f or simplicity we disregard intervalley 
scattering in our discussion of optical phonons below. 
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IV. ACOUSTIC PHONONS IN MONOLAYER 
AND BILAYER GRAPHENE 

Long-wavelength acoustic phonons may be treated 
with continuum theories to a good approximation. An 
atomistic treatment is required only for optical phonons, 
see Sec. lVIl Here we do not discuss the actual calculation 
of the phonon modes, which could be described in terms 
of elasticity theoryP^^U 



A. Coupling via deformation potential 

The dominant form of electron-phonon coupling for 
long-waveleng th acou stic modes is due to the deforma- 
tion potential! 17 * 25 * 27 * For in-plane modes the potential is 
of the form DV u, where u(r) is the displacement vector 
from elasticity theory, and D is the coupling constant, for 
which values in the range D = 10 — 50 eV have been used 
(Refe.EEE9. 

The deformation potential is nonzero only for longitu- 
dinal (LA) modes. Thus we may use the expansion in 
terms of longitudinal plane waves: u(r) = J2 q u q qe lq r . 

Quantizing this with u q = i^/H 2 / (2Mu q )(b q + b_„), 
where b q and b q are the phonon creation and annihilation 
operators, the electron-phonon coupling Hamiltonian be- 
comes 



right-hand side describes a conical surface in momentum- 
energy space with the apex at (fc, e£) and the left-hand 
side is another surface. In the case of MLG this surface is 
also a cone, while for BLG in the parabolic-band approx- 
imation it is a paraboloid. The allowed p values now lie 
on the p-planc projection of the curve defined by the in- 
tersection of the two surfaces. These may be worked out 
analytically, and they are needed for the numerical eval- 
uation of the heat current. The details are given in App. 
|A| It is noteworthy that in the case of MLG there cannot 
be any interband scattering (/3 ^ a) when the phonon 
dispersion is linear,^ while for BLG interband scatter- 
ing is possible if the initial k is small enough. However, 
in practice this scattering is strongly suppressed by the 
coupling constant (13). 



V. HEAT CURRENT BETWEEN ELECTRONS 
AND ACOUSTIC PHONONS 

Here we use the results of Sec. [H] for graphene, with 



the LA phonon coupling constant given by Eq. ( 13 ). First 
we present the general result which can be evaluated nu- 
merically, and then discuss analytic approximations in 
the BG and EP limits. 



Numerical solution 



ka,p@ q 



(12) For acoustic LA phonons the most general expression 

for Qi n d without any approximations is 



with M^ q = -^/h 2 /(2Muj q )Dq(pP\e iq - r \ka). Here 
Lu q = hcq is the LA phonon dispersion and M = Ap is the 
total mass, with p the mass density (of MLG or BLG) 
and A the area of the system. The wave functions are 



(r\ka) — e ' r 4>%/V A, where Eqs. (11) should be used 
for The coupling constants are then identified from 



«Op.*+«' which Yield**™ 



w 



a/3 = n 

k P'i ~2Moj, 



D 2 q 2 F aP (8). 



(13) 



Here F a p(Q) = (1 + s a sp cos n0)/2, 6 = <j) p — and 
n = 1,2 for MLG and BLG, respectively. This shows 
that for MLG backscattering from acoustic phonons is 
weak, whereas for BLG there is no difference between 
forward and backward (intraband) scattering. 



B. Conservation laws 

In scattering of electrons from acoustic phonons, the 
final wave vectors p allowed by energy and momentum 
conservation laws (for given initial fc) may be determined 
from the equation = e£ + shc\p — fe|, where s = +1 
for absorption and s = — 1 for emission, and and 
are the initial and final energies, respectively. The 



,AhD 2 



E 



1 + s a s y(k,q) 



OO /• <?max 

dkk I dqq 2 



^l-{x{k,q)) 2 ' 
xn ac (ncq)[m-hcq)-f(e%)], 



(14) 

i-2|-l v ' 



with n = 1,2 as above. Here v% = s a n(hv) n (k/ji) n ~ 1 jh 
is the Fermi speed and x(k,q) = — {[^(fc,^)] 2 — k 2 — 
q 2 }/(2kq) is the cosine of the angle between the in- 
coming k and the phonon q. In the latter p(k, q) = 
{s a sp — sphcq r yi~~ 1 /(hvk) n ) 1 / n k is the length of p = k — q 
after imposing conservation laws. Moreover, y(k, q) = 
n(z{k, q)) n — n+l, where z(k, q) = [k — qx(k, q)]/p{k, q) is 
the cosine of the angle between fc and p. Finally, the cor- 
rect limits qZ!t v of the q integral are obtained from the 
conservation laws, as discussed in Sec. |IVB| Using these 
results together with Eq. the total power Q e -ac may 
be obtained numerically. An example of the correspond- 
ing thermal conductance G e _ ac (/^,T) is represented by 
the solid lines in Fig. [2j 



To obtain analytic estimates, we first expand Eq. (14) 
to leading order in c/v c 1. In this approximation we 
may neglect interband transitions (j3 ^ a) also for BLG 
and use x{k,q) — q/2k, p(k,q) = k, q^ K = 2k, and 
q^? n = 0. In the BG and EP limits the forms of Eqs. 
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FIG. 2: (Color online) Thermal conductance 

G e - ac (n,T) for MLG (a) and BLG (b). The unit 
is Go = g e AD 2 h\ f i\ 4 k B /[8n 2 p{hv) (i ] for MLG and 
Go = gcAD 2 hj 1 \fi\ 3 k B /[l6TY 2 p(hvf] for BLG. The pa- 
rameters are c/v = 0.02 and *ya / 1 A* I = 10 m 3 )- The solid 
lines (magenta) indicate the full numerical solutions, and the 
dashed lines (black) the analytical approximations. 



tronic specific heat is given by C — (ir 2 /?>)v(p)k 2 B T ', 
where the density of states including the degeneracy 
g e = 4 is u(E) = 2\E\/{nh 2 v 2 ). Identifying r^" 1 = a*T 2 
and C = 7*T, we find £ = (2/5)a* 7 *. 

Seeking a high-temperature approximation, we re- 
quire hcq <C k B T e ,k B T ac , which allows us to expand 
also the Bose function. Since for a given k the max- 
imal q is 2k, and the maximal relevant k is deter- 
mined by \e^\ ~ max(|/i|, k B T e ), the limits translate to 
2(c/v)max(\p\, k B T e ) -C kBT e ,kBT ac . This is the case 
considered in Ref. El and since c/v <C 1 the limit is equiv- 
alent to T e ,T ac > T bg ,mlg, if T e ~ T ac . In this limit, 
Qind depends on both T e and T ac . Using Eq. Q then 
leads to a total power of the asymmetric form in Eq. ([2| , 
where Vd = A and 



D 2 k B 
~30irph 5 v 6 
+ 7n 4 (k B T [ 



15p 4 + 30ir 2 p 2 (k B T e ) 2 



(17) 



([!]) and ([2| are then recovered. Below we give the results 
for the coefficients £ and g(fj,,T e ) for MLG and BLG. 
The thermal conductances G e -ac obtained from them 
are shown as the dashed lines in Fig. [2| The agreement 
with the full numerical solution is very good. The most 
notable difference between MLG and BLG is the much 
weaker temperature dependence of G e _ ac for k B T e 3> 
in BLG. 



B. Limits for monolayer graphene 



For fcB?BG,MLG ^ k sT e < \p\ this yields g(fX,T e ) oc p 4 , 
and in the limit \p\ <C k B T e we find g(fi,T e ) oc T 4 

For completeness we mention also the result!^ for the 
high-temperature relaxation rate in MLG: 



D 2 \p\k B 
2ph?v 2 c 2 



In this case the relation to Q e 
than at low temperature. 



T. (18) 
is more complicated 



First consider the case of MLG (n = 1). For a low- 
temperature approximation of Eq. ( |14[ ) we may estimate 
£fc = s a hvk — (j, and hcq — k B T ac . If k B T e ,k B T ac < 
2(c/v)\p\ = fcs Jbg,MLGj w e may set q/2k — >• in the 
g integral and extend the upper limit q max = 2k — > oo. 
Since Ticg -C |e^|, we expand the Fermi function to the 
first order in hcq. The T e dependence then drops out of 
Oind so that the total power takes on the symmetric form 
of Eq. ([I]), where = A, S = 4, and 

15pfi 5 i; 3 c 3 ' V j 

The same result was derived previously in Ref. [H 

The coupling constant £ may be expressed in terms 
of the electron-phonon relaxation rate at the Fermi level 
and the electronic specific heat P We start from the defi- 
nition of the relaxation rate t^ 4 = —df%S e -ph{fk)\f=fa> 
where f (E) = {exp[(£ - p)/k B T] + l}" 1 and T is the 
lattice temperature. A similar calculation as for Qi n d 
yields 



nD 2 k 2 9 



Aph?v(? 



(16) 



where hp = \p,\/(Hv). (We n ote t hat unlike this result, 
the transport relaxation rat d 14 * 17 l is oc T 4 .) The elec- 



C. Limits for bilayer graphene 

Next consider BLG (n = 2). In this case, besides p, T e , 
and T ac , there is an additional energy scale determined 
by 71 . In order for the parabolic two-band approximation 
to be valid, we must require max(|/z|, k B T e ) -C 7l . 

To find a low-temperature approximation, we follow 
similar steps as for MLG. The assumed limit is now 
k B T e ,k B T ac < 2(c/v)^j 1 \p\ = fc B T BG ,BLG (with |/z| < 
7i), using which we find Eq. ( T]) with 5 = 4, but this time 



S = 



v 2 D 2 lik 4 B 
60ph 5 v 3 c 3 



(19) 



The scaling with temperature is thus similar to the MLG 
case, but the p dependence is different. Result (19) is 



valid also in a normal 2-D system with effective mass 
ll /(2v 2 ). 

For BLG we find the relaxation rate 



irD 2 k 2 B 
8ph 3 vc 3 



(20) 



where kp = \J^j\\p\/ {hv). (The corresponding transport 
relaxation rate is again oc T 4 .) Now the density of states 
entering C is approximately v(E) « 71/ (irh 2 v 2 ). Also in 
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this case the definitions r k 1 = a*T 2 and C — j*T lead 
to S = (2/5)a*7*. 

For BLG the high-temperature (EP) limit 
hcq <C k B T ei k B T ac may be written as 
2(c/w) v / 7imax(|^|, k B T e ) < k B T e ,k B T ac . This is 
again equivalent to T e ,T ac ^> T B g,blg, if additionally 
T e ~ T ac ^ 2(c/w) 2 7i/fcB- The total power again takes 
the form of Eq. pi , where now 



D'ljk B 
4irph 5 v 6 



{2k B T e \n[2cosh(p/2k B T e )}} . (21) 



In contrast to MLG, if /cbTbg.blg "C fcsT e <C then 
g(n,T e ) (x p, and if |^|, 2(c/v)^jik B T e < fc B T e then 
g(p,T e ) oc fc B T e . 

Finally we mention also the BLG result for the high- 
temperature relaxation rate, 



D 2 ll k B 



(22) 



k " 4ph 3 v 2 c 2 ' 
which is valid in the same limit as assumed above. 

VI. OPTICAL PHONONS 

For graphene, its multilayers, and graphite the phonon 
spectra have been studied in detail both experimentally 
and theoreticallyP^SI i n order to describe the crossover 
from acoustic phonons to optical phonons as the domi- 
nating phonon type for heat dissipation, we use some sim- 
plified optical-phonon models. For the description of the 
intrinsic optical modes an atomistic description is needed 
as a starting point. However, we skip the details here 
(see App. [B]), as similar calculations have been reported 
earlier. Of the intrinsic phonons, we consider explicitly 
only in-plane LO and TO modes (collecti vely LT) at long 
wavelengths, i.e. near the T pointP^^Eim! (The K-point 
intrinsic phonons^ are expected to have coupling con- 
stants that differ only by a numerical prefactorPS) Addi- 
tionally, we consider co upling t o the "remote" phonons 
of a dielectric substrate) 15 * 37 * 38 * which is also a relevant 
concern for many experiments.^ 



A. Coupling constants for simple phonon models 

Let us first consider a model for the long-wavelength 
in-plane (LT) optical phonons in MLG, for which the 
LO and TO branches are nearly degenerate, with energy 
Ult ~ 0.2 eV. In the tight-binding picture of Fig.[T] the 
strongest coupling to the LT modes comes from the mod- 
ulation of the nearest-neighbor (Al-Bl) coupling ampli- 
tude _7o ( see Ref- EH)- The coupling constant is of the 
forrrPl 



w 



iLT 9(%) 2 h 2 1 



hp.q 



2MCI 



LT 



(l-s a S0 cos(0 fc + </>p-20a)), (23) 



where M — Ap%, with p\ w 7.6 • 1CP 7 kg/m 2 the mass 
density of MLG, and 7q ~ 40 eV/nm is the derivative of 
7o with respect to the nearest-neighbor bond length.^ 
The vector a(q) = q for the LO and a(q) = z x q for the 
TO branch, z being normal to the plane. 

In BLG there are four nearly degenerate LT branches, 
with Qlt ~ 0.2 eV, since for both LO and TO type 
modes the atoms in layers 1 and 2 can move either in 
phase or in opposite phasesPf With perturbation theory 
to first order in hvk/ji we find in this case 

a ,, LT _ 9(%) 2 h 2 l(hv) 2 ( 

w k P , q - 2 Mn LT 2^~\ k +p 

+ 2kp [± cos(0 fcp ) - s a s p cos(0 fc a + P a)] 
T s a Sj3 [k 2 cos(2<f> pa ,) + p 2 cos(20 fed )] 

' (24) 

where M = Ap2, P2 = 2pi, and 4>kp = 4>k — 4>p- The 
upper and lower signs are for the in-phase and opposite- 
phase modes, respectively. For BLG there are also the ZO 
modes with flzo ~ 0.1 eV which could couple linearly to 
electrons via the modulation of the A1-B2 bond and thus 
7i , but we estimate them to be relatively unimportant for 
the energy relaxation (see Ref. |4"0"|). 

Finally, in addition to these intrinsic phonons, we con- 
sider possibility of co upling t o remote surface phonons of 
a dielectric substrate! 15 * 16 * 38 ! The coupling is due to the 
electric polarization associated with the phonons, which 
modulates the scalar potential on the graphene. For sim- 
plicity we again disregard the possibility of intervalley 
scattering.^ The coupling constant is given b} * 15 * 37 ! 



a/3, rem 
^kp,q 



e Ml T 



1 



2e A q 



- 2qz F aP {6), (25) 



where £l rcm is the energy of the relevant surface mode, 
z > is the effective distance between the graphene and 
the substrate, F a p(9) is as in Eq. (13), and £o is the per- 
mittivity of vacuum in SI units. If there is only a single 
relevant surface mode, /3 rom = (e s — £ tx) )/[(e s +l)(e 0O +l)], 
where e s and are the static and high-frequency di- 
electric constants of the insulator.^ If there are several 



surface modes, Eq. (251 should be additionally weighted 
by the corresponding relative oscillator strengthsP^ 



B. Heat current between electrons and optical 
phonons 

When the momentum conservation is imposed, q = 
\Jk 2 + p 2 — 2kp cos 6, where 9 = <pk — <i>p- All the cou- 
pling constants may then be expressed in terms of only 
k, p, and cos 9. For an electron dispersion of the form 
€ k = s a£k, £fc > 0, the heat current between electrons 
and dispersionless optical phonons may be written 



)(7) 
* e—op 



Aq^J (/i,T e )[n e (fi 7 ) - n ph (fl 7 )} 



(26) 
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where (including degeneracy of g e = 4) 



de-op m 



=— Vtf 



f(r2 7 x)z/(fi 7 (a; — 1)) 



x / d9w al3 ^(k(QJx\),k(n^x- 1|), 



COS ( 



x [/(Q 7 (z - 1)) - f{^x)]dx. 



(27) 



Here a/3 are chosen such that s a = sign(cc) and sp = 
sign(x — 1), f(e) is the density of electronic states, and 
k(e) the inverse of efc. The interval < x < 1 corresponds 
to interband scattering. 

In the case of the LT modes further simplification is 
achieved by summing over the degenerate LT modes 7. 
The angle- dependent terms of the coupling constants ( 23 ) 
or (24) then cancel and the angle integral in Eq. (27) 



becomes trivial. Thus for MLG' 



ir(hv) 4 pi 



(28) 



where the factor g e — 4 is included, SI = S1_lt, and 

/oo 
\x(l-x)\[f(ft(x-l))-f(nx)]dx. (29) 
-oo 

For \i — and f2/^B?e » 1 it has the value T = 1/6. 
Similarly for BLG 



{LT) mvH^fhj, 

ql_ ^,T e ) = — — y —G (fi,T e ), 



(30) 



where 



S(/x,T e )= / -(|x| + |x-l|)[/(^-l))-/(^)]dx. 

J —oo 

(31) 

For = and Sl/fcsT e ^ 1 this function has the value 
Q = 1/4. We note that due to the assumed parabolic 
dispersion, this result is a good approximation only close 
to = 0. For the case of the remote phonons, the angle 
integral in Eq. (27) is quite complicated, but in the limit 
2 — it may be carried out analytically for both MLG 
and BLG. However, the expressions are complicated and 
we skip them here. 



The numerical results for the prefactors q. 



(t) 



-op 



arc 



shown in Fig. [3] For the remote phonons we use z — 
and the values for 6H-SiC (see Ref . |4"T|) . where e s = 9.72, 
£oo = 6.52, il rcnl = 116 meV, and hence /3 rom = 0.040. 
It is seen from Fig. [3] that the prefactors for the remote 
phonons are considerably higher than those of the in- 
trinsic phonons. This fact combined with fl lcm < Qlt 
means that for heat dissipation the remote phonons may 
in practice always be the more important ones. The limit 
z = overestimates the coupling constant somewhat, 
but 9e-op decays with distance at most proportionally 
to e _29maxZ , where for example for MLG at k-oT e <C Orem 




0.05 



S LT ) 



and g^!!™ of the 



FIG. 3: (Color online) Prefactors q e _, ,, ,,,, 

electron-optical phonon heat currents for different doping lev- 
els, (a) and (c) are for MLG, and (b) and (d) for BLG. Upper 
panels are for the intrinsic LT phonons and the lower ones for 
the remote phonons, with 2 = 0. 



we find g max = (Sl rcm + 2\/j,\)/(Hv) and thus q~l x rj 6.2 
nm when /i = 0. The use of Si02 would reduce e s — e^, 
but then one of the two modes with appreciable oscilla- 
tor strengths also has a lower energyPS The correspond- 
ing crossover temperatures between Q e -ac and Q e - op are 
discussed below. 



VII. DISCUSSION AND CONCLUSIONS 

Let us briefly discuss the implications of our results. 
A common experimental situation where the heating of 
electrons occurs is that of a two-probe measurement with 
a finite bias voltage V. The Joule heat created at some 
point of the system is dissipated from the electrons in ba- 
sically two ways, either by a transfer to phonons or by dif- 
fusion of the electrons away from the heated region (e.g., 
into the electrodes). The static situation is described by 
a heat balance equation V • jg + -P c -ac = Pj, where jg 
is the heat current density, while P - ac and Pj are the 
local electron-acoustic phonon power and Joule power 
per area, respectively. We assume for simplicity that 
the Wiedemann-Franz law applies (see Ref. 142]) so that 
Jq = — «VT e , where k — CaT e is the heat conductiv- 
ity, a the electrical conductivity and C — (7r 2 /3)(/cs/e) 2 
is the Lorenz number. To find some order-of-magnitude 
estimates, we assume a quasi-one-dimensional situation 
with a sample of length L, and consider the simpler, dis- 
cretized equation P^iff + Pc-ac = Pj with the diffusion 
power P diS = (4£a/L 2 )(T 2 - T 2 ) and Joule power Pj = 
aV 2 /L 2 . We also concentrate on the low-temperature 
regime, k B T < fcs?BG,MLG = 2(c/v)| / u| for MLG and 
k B T < /j b Tbg,blg = 2(c/u)V7iR for BLG, so that 
P c - ac is of the form P c _ ac = E(T 4 - T 4 ). Here T is the 
electron temperature in the middle of the graphene sam- 
ple and we assume the acoustic phonons and the electrons 
in the leads to remain at the bath temperature To- 

Clearly, at low enough Tq and bias V the diffusion 



power dominates over the phonon power (P c - ac < -Pdiff)- 
However, with increasing To or V there is a crossover to 
an electron-phonon dominated regime (P c _ ac > Pdiff). 
If T > T 0x — 2£cr/(SL 2 ) then P e _ ac dominates also 
in the linear-response regime (arbitrarily small V and 
AT = T — To). For To < To x there is a finite crossover 
voltage V a = \/8£(T c 2 r — Tq ), with the corresponding 
temperature T cr = y/4Ccr/(T,L 2 ) — T 2 . 

When the phonon power dominates, P -ac ~ Pj- Since 
Pj can be deduced from the current-voltage character- 
istics, measurement of the electron temperature in the 
presence of heating can also act as an indirect measure- 
ment of the electron-phonon coupling constant. To see 
if this regime can be reached before the low-temperature 
approximation for P c -ac breaks down, we estimate the 
crossover temperatures. Using D = 30 eV we find for 
MLG 



a 



a 





FIG. 4: (Color online) The ratio of (LT) optical phonon ther- 
mal conductance G'^^/i, T) and that of acoustic phonons 
G e -a C (p,T) for MLG (a) and BLG (b). The parameters are 
D = 30 eV and \ j' \ = 40 eV/nm (Ref. H7J|. In both cases the 
crossover point G^jL/Ge—ac = 1 moves to higher tempera- 
ture with increasing 



and for BLG 



5 K x 



Woo) 



1/2 



T 0cc = 11 K x 



(L/l A im)(|/i|/0.3 eV)Va 



(a/ao) 1 / 2 (| M |/0.3 eV) 1 / 4 
(L/l M m) 



(32) 



(33) 



where uq = 4e 2 /h. These are only valid if Tq x <C 
Tbg.mlg/blg- We find T bg ,mlg = 140 K x (M/0.3 eV) 
for MLG and Tbg,blg = 160 Kx (M/0.3 eV) 1 / 2 for BLG. 
Since the electron-phonon interaction becomes more im- 
portant for increasing T, it seems possible to meet these 
criteria with long enough samples. 

The estimate could be improved by taking into ac- 
count that a also grows with /1 differently for MLG and 
BLG, and it may also depend on temperature! 43 * 44 ^ How- 
ever, these details depend on the types of scattering. In 
practice, the dominant form of scatt ering in graphene 
on a substrate is from impurities ! 45 * 46 ! Experiments indi- 
cate that a grows roughly linearly with charge density,^ 
which yields a ~ /i 2 for MLG and a ~ /x for BLG. We 
also note that the diffusion power can be reduced or even 
eliminated by the use of superconducting leads. 

At high enough temperature there is another crossover, 
where optical phonons beg in to dominate the electron- 
phonon heat transfer P1°J Although the conductance 
G e _ op is exponentially suppressed at low tempera- 
ture, G e -ac is also small and for MLG quite strongly 
temperature-dependent, with G e _ ac ~ T 4 at /1 = 0. 
This can lead to surprisingly low crossover temperatures. 
Here we consider explicitly only the intrinsic LT optical 
phonon modes. Figure 4 shows the ratios G^^ v jG e - ac 
of optical and acoustic phonon ther mal c onductances for 
MLG and BLG. As found previouslyp^ optical phonons 
can become dominant already well below room temper- 
ature (150 K for MLG and 200 K for BLG), and the 
crossover moves to higher temperatures with increas- 
ing carrier density. The non-monotonous behavioral of 

( LT) 

G e _'/G e - ac in the case of MLG is not present in BLG. 



As found above, if the graphene is on a dielectric sub- 
strate, the dominant optical modes are most likely the 
surface modes of the dielectric! 15 * 16 ! By using the param- 
eter values quoted in Sec. |VIB| for SiC and studying 
Gg r ™p /G e -ac as in Fig. 4j we find crossover tempera- 
tures 60 K for MLG and^O K for BLG at fi = 0. For 
Si02, using the parameters quoted in Ref. [TBI the results 
are as low as 30 K and 50 K, respectively. 

We note that the above results depend on the poorly- 
known parameters D and 7 in their second powers, while 
c and v also appear with high powers. This makes quan- 
titative predictions difficult, and increases the need for 
an experimental determination of the coupling constants. 
Finally it should also be noted that if the bias voltage 
exceeds the value ~ Oir/e = 0.2 V, scattering from op- 
tical phonon s becom es very important also at low bath 
temperature j 16 I 48 ]H1 j n this case the optical phonons can 
have a highly non-equilibrium distribution! 5 -^" 5 - 2 ^ 

To conclude, we have calculated the power transfer be- 
tween the electron and phonon systems in monolayer and 
bilayer graphene, assuming the existence of quasiequilib- 
rium. In particular we have studied the coupling to lon- 
gitudinal acoustic phonons and different types of optical 
phonons. For the former we have calculated the power 
numerically and derived analytic expressions in low- and 
high-temperature limits. The power transfer to acoustic 
phonons dominates the diffusion power above a certain 
crossover temperature, estimated in Eqs. (32 1 and (33 1. 



At even higher temperatures, there is another crossover 
where optical phonons begin to dominate, and we have 
estimated also these crossover temperatures numerically 
(Fig. [4]). We find that for graphene on the substrate the 
most relevant optical phonons are likely to be the surface 
optical modes of the substrate. 
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FIG. 5: Momentum conservation in absorption (+) or emis- 
sion (— ) of a phonon. fc is the incoming electron wave vector, 
p± the outgoing one, and q that of the associated phonon. 




FIG. 6: (Color online) Schematic illustration (c/v not to 
scale) of the allowed final states in the momentum-energy 
space for electron-acoustic phonon interaction in the case of 
MLG (a) and BLG (b). These are given by the curves along 
which the two momentum-energy surfaces intersect (thick 
solid lines). The curves corresponding to elastic scattering 
are also shown (thin solid line). The initial state (with fc 
in the direction of the arrow) is in the conduction band and 
only the case of emission is shown. In the case of MLG, no 
interband scattering is possible. 
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FIG. 7: (Color online) Allowed final states in the p plane for 
electron-acoustic phonon interaction in the case of MLG (a,b) 
and BLG (c,d) graphene. The upper diagrams (a,c) show the 
intersections of Fig. [6] in (p\\,E) plane, where we define py, 
p± as the components of p parallel and perpendicular to fc. 
The lower diagrams (b,d) show the curves in p plane. The 
thin lines correspond to elastic scattering, thick solid lines 
to emission, and thick dashed lines to absorption. Interband 
scattering is possible only for BLG, and here only for emission 
(initial state in the conduction band). For MLG we use the 
unphysically large c/v = 0.05 for clarity of the figure, and for 
BLG x c {k) = x c (k cl ) + 0.005 = \/2 - 1 + 0.005. Energies are 
in units of Ek = hvk. 

curves are then described parametrically by the functions 
q — q(k, cos </>). We give the solutions first for MLG and 
then for BLG. The results are illustrated in Fig. [7] 



Appendix A: Analytic solutions for the conservation 
laws in scattering of electrons from acoustic phonons 

Here we detail the solutions for the possible final states 
allowed by momentum and energy conservation laws in 
electron-acoustic phonon scattering. If k and p are the 
initial and final wave vector of the electron, respectively, 
q the wave vector of the associated phonon, and k = \k\, 
p = \p\, and q — \q\, then the conservation laws are 



P 

c . 



k 



sq, 

- shcq. 



(Al) 



Here s = +1 for absorption and s = — 1 for emission, 
and a = c, v is the band index. The final wave vec- 
tors p s = k + sq are illustrated in Fig. [5] Figure [6] il- 
lustrates schematically the allowed states in momentum- 
energy space. 

We represent the final states in polar coordinates by 
writing p s = k + s(q cos <j)p\\ +(?sin0p ± ), where py is the 
unit vector parallel to k and the one perpendicular to 
it, and <f> is the angle between q and k (see Fig. [5]). The 



1. Monolayer graphene 

For MLG e£ = s a hvk, with s a — ±1. Eliminating p 
in Eq. (Al), we find 



sptivp s (k, q, x) = s a hvk + shcq, 



(A2) 



where p s {k,q 1 x) — \p s \ — \Jk 2 + q 2 + 2kqsx, with p s = 
k + sq (s = ±1), x = cos0 = q ■ k/(kq), and s = ±1 
for absorption or emission, respectively. This yields the 
parametric representation for the curves 



— sx+ss. 



.(c/v) 



q_ 

2k 



i- {c / v )2 , sx < ss a (c/v), P = a 
0, sx > ss a (c/v), P = a (A3) 

no real solution, j3 ^ a. 



There are thus no finite-q solutions if sa; > c/v, one 
for — c/v < sx < c/v and two if sx < —c/v, so that 



lim. 



c/v. The maximal values are (set sx = — 1) 



9max _ I l-ss a (c/v) ' (3 — a 



2k 



no real solution, /? 7^ a, 



(A4) 
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while the minimal value is always q^ n — 0. In the 
special case c/v = the scattering is elastic and q = 
—2ksx, which has the maximal value 2k corresponding 
to backscattcring. If s a = +1, then q^ x < 2k for emis- 
sion (s = —1) and q^ x > 2k for absorption (s = +1), 
and vice versa for s a = —1. 



2. Bilayer graphene 

For a BLG e£ = s Q (/wfc) 2 /7i, with s a = ±1. Now 

s^hvp^k.q.x)) 2 = Saihvk) 2 /^ + shcq, (A5) 



where again p s (k,q 7 x) = \Jk 2 + q 2 + 2kqsx, with x = 
cos0 = k ■ q/(kq). Now we find the following solutions. 
For f3 = a 

—sx + ss a x c (k), sx < ss a x c (k), k > k C 2 
q I 0, sx > ss a x c (k), k > k c2 

2k I — sx + x c (k), < k < k C 2- ss a = +1 
0, < k < k c2 , ss a = -1, 

(A6) 

where x c (k) = (l/2)(c/v)[^i/hvk) and k C 2 = 
(l/2)(c/w)(7i//iw), so that x c (k C 2) = 1. For k > k c2 
there are thus no fmite-g solutions if sx > x c (k), one for 
x c (k) < sx < x c (k) and two if sx < — x c {k), such that 
liniq-^o |a;| = x c {k). No solutions exist for k < k c ±. The 
extremal values of q are 



/3=a 
?max,min 

2k 



1 + ss a x c (k), k > k c2 
±1 + x c (k), < k < k C 2, ss a = +1 
0, < k < k C 2, ss a = — 1, 

(A7) 



where the maximum (minimum) follows by setting sa; 



T (sa; 



+1). Again, q. 



0=a 



2k for absorp- 



tion/emission, if s a — +1 and vice versa for s c 
For (3 ^ a we have 



•1. 



q_ 

2k 



-sx + x c (k) ± \j (—sx + x c (k)) 2 



(A8) 



if |sa; — x c (k)\ > y/2, < k < k c i, ss a = —1 where 
k c i = k C 2/(v2 — 1) > k c2 . For k > k c \ or ss a = +1 no 
real solution exists. In this case there are only solutions 
for either emission (s = —1) or absorption (s = +1), 
depending on s a — — s^ such tha t ss p = — ss a = +1. 
The upper and lower signs in Eq. ( A8 1 describe the two 



branches of the interband solution, represented by the 
closed circular curve in Fig. [7} By setting sa; = —1, we 
find that the extremal values for q are now 



?max,min 

2k 



1 + x c (k) ± ^(l + x c (k)) 2 -2 



(A9) 



if < k < k c \ and ss Q = —1. We also note that 
x c (k cl ) = V2 - 1 w 0.41421 and thus q^ a j2k cl = 




FIG. 8: (Color online) Direct and reciprocal lattices for 
graphene, showing the primitive vectors ai_2, bi,2 of the Bra- 
vais lattices and the nearest-neighbor vectors ^1,2,3. For the 
reciprocal space the high-symmetry points F, K, K', and M 
are shown, and the hexagons around the T points are the 
first Brillouin zones. The two inequivalent parts of the Fermi 
surface are sketched with circles. 



9clfmin/2fccl 

value for the intraband solution (ssp 
9mlx. c i/2fcci = 2 - V2 « 0.58579. 



1/V2 w 0.70711. The corresponding 

ss„ = —1) is 



Appendix B: Atomistic description of 
electron-phonon coupling in graphene 

Here we consider the tight-binding description of the 
coupling of electrons to the intrinsic phonons in graphene, 
and explain in more detail our simple models for the long- 
wavelength optical phonons in MLG or BLG. We only 
consider BLG in detail. 

For the unperturbed system (without electron-phonon 
coupling) the tight-binding Hamiltonian of BLG is of the 
form 



(Bl) 



where = (a\ k , b\ k , b' 2k , a\ k ) consists of the electron 
creation operators at different sites and 



H(k) = hvk 



( 

ji/Hvk 

\ 



±e ±i* 







^i/ hvk 



_|_g±i0fc 










\ 



(B2) 



with ipk = arctan(fc y /A; x ) and hv — 3a7o/2. The upper 
and lower signs are for the K and K points, respectively 
(see Fig. [8|. The zero of energy is thus set at the charge 
neutrality point. The matrix H(k) has four eigenvalues 
(energy bands) and eigenvectors, which we den ote e a (k) 
and $ Q (fe), a — 1,2,3,4. The Hamiltonian (Bl) and 



other similar operators may be written in the eigenba- 
sis by expanding ^j, = J2 a ^a(k)ck a where cj~ a is the 
annihilation operator for the eigenstates. 
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1. Electron-phonon coupling Hamiltonian 



we find 



The coupling to the in-plane modes (here LO,TO) is 
obtained by assuming that when the equilibrium nearest- 
neighbor distance a is perturbed by 6a, then 70 — > 70 + 
7o<Sa. For BLG this lead s to an electron-phonon coupling 
Hamiltonian of the for d 12 ' 34 ' 36 ' 



^-S=tE E 1&i(H)lM*+*i) 

K 3=1,2,3 
X 8 3 ■ [UAl(R) -U B l(R + Sj)] 

+ ^E E ^A2(R.)i>m{R. + Sj) 

R s 3=1,2,3 

x 5j • [w^ 2 (-R s ) - u B 2(R s + Sj)] + h.c, 

where R sums over the Al lattice sites, R s are the A2 
sites, and where u$ a is the in-plane displacement of atom 
5 = A, B in layer a = 1,2. Here we may insert the Fourier 
transformations 

i) A a(R) = -j= e <fe "Vfe, <r = l,2, (B4a) 
^ fc 

w Aff (H) = i E eU, ' Hu A*v * = !. 2 » ( B4b ) 



(B3) 



where iV is the number of unit cells, and similar ones 
with ^4 — > B and a CT fc — » fro-fc. Then, calculating the wave 
vectors relative to K or K' and assuming q <C = 
47r/(3v / 3a), we obtairP 



i? 



37o 1 



(LT) _ 

e—ph 2 








EE*l 



p+<r 



U y lq T iUf q 





UZ,TiU% 






U|g±^2 X g 



2q 



(B5) 



where 17^ = U Aaq — U Bag (f = 1,2) is the relative 
displacement vector between A and B atoms. 

For the ZO optical modes, which can modulate the 
equilibrium distance d between the Al and B2 atoms 
such that 71 — > 71 + j[8d, the coupling Hamiltonian is of 
the form 



H 



(zo) _7i 



e—ph 



d 



J2 E i^Al(R)i>B2(R d ) 
R 3=1,2,3 

VLW^Al(i*)PB 2 (-R<i) ~ /lAl(fl)] 



(B6) 

Here hs a is the out-of-plane displacement of atom 6a and 
i?d is the site of the B2 atom above R. Inserting Eq. 



(B4a) and 



aq 



(B7) 



(ZO) 
e—ph 



7i 



/ ffg 0\ 



i? q 

V 0/ 



P" 



(B8) 



where = H B 2 q - H A i q . 

These results must still be expressed in the low-energy 
eigenbasis by using the expansion "J/- = J2 a ^a{k)ck a , 
and where a sums only over the low-energy bands (c, v), 
keeping only terms to first order in Hvk/ji. Also, a trans- 
formation of the form 
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(B9) 



where C — x iVi has to be introduced, with appro- 
priate models for the matrix X(q) as discussed be- 
low. Finally, the phonons are quantized with U ql — 
y /h 2 /(2il 1 (q))(b q7 + &Lq 7 ), where fi 7 (g) is the energy of 
the phonon mode. This leads to the form 

He-ph = E E E M pk'2 C l^{K, + &-g-y) (B10) 
Q7 ka p/3 

such that the matrix elements M^' 7 and thus the cou- 
pling constants may be identified. For MLG the deriva- 
tion is similar. 



2. Coupling constants for simple intrinsic phonon 
models 

Next we consider some simple models for long- 
wavelength optical phonons! 27 * 34 * 35 ^ We do this by 
specifying the form of the transformation matrix 
X(q) for the different phonon branches, which should 
otherwise be obtained by diagonalizing the phonon 
Hamiltonian. The matrix must be normalized as 
£<5aC X *SaC, 4 X ( l) M s<y x SaC,M = 6^, where M &a is the mass 
of atom 8a in the unit cell. In graphene M$ a = Mc, the 
mass of a carbon atom. Below we use the vector notation 

Let us first consider a model for the long-wavelength 
in-plane (LT) optical phonons in MLG, for which the 
LO and TO branches are nearly degenerate. To describe 
the characteristic opposite-phase motion of the A and B 
sublattices we choose 

x ^ lt(<7) = vm &(q) ' XmLT{q) = vm 6{q) ' 

(Bll) 

where a(q) = q for the LO and a(q) = z x q for the TO 
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branch, z being normal to the plane. This yield^^l 

w kp, q = ^iX-s a spcos[(j) k + ^ p -2<t> & )), 

(B12) 

where M = 2M C N. 

In BLG there are four nearly-degenerate LT branches. 
For them we choose similarly 

X A1 , LT (q) = —=a{q), X m>LT (q) = —==a{q), 
zy J-Vlc iy Niq 



Xa2,lt{q) = — 7 ==a(q), 
£\J Mc 



X b2,lt{q) = -^—a{q). 

(B13) 



According to these, for both LO and TO type modes the 
atoms in layers 1 and 2 can move either in phase (upper 
signs) or in opposite phases (lower signs). We find 



kp,q 2M9. L T 2 7l 2 

+2kp [± cos(0 fcp ) - s a s/3 cos(0 fe a + P a)] ( B14 ) 

Ts a Sfi [k 2 cos(2(j) p& ) + p 2 cos(20 fcA )1 



with M = 4M C N and 4> kp = <j>k~<f>p- To model the ZO 
modes of BLG, we choose 



Xai,zo(q) 
Xa2,zo(q) 



1 



2 V / Mc 
1 



z, x B \,zo{q) 



-l 



2VM 



B2,ZO 



(9) 



2VMc 
-1 



2VM 



(B15) 



where the signs amplitudes of the A2 and Bl displace- 
ments are unimportant in our approximation, which ne- 
glects any coupling between them. This yields the cou- 
pling constants 



a p,zo 2(7^) 2 S 2 2 



(B16) 



with M — AMcN. This form is particularly simple, be- 
cause there is no angular dependence, and no interband 
coupling. In fact, for k — p (q — >• 0) the result may be 
obtained from the continuum description by modulating 
7l,ie. the effective electron mass 7i/(2v 2 ). 
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